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We consider the proximity effect in a Topological Insulator sandwiched between two conven- 
tional superconductors, by comparing s — wave spin singlet superconducting pairing correlations 
and odd — parity triplet pairing correlations with zero spin component orthogonal to the slab ("po- 
lar " phase). A superconducting gap opens in the Dirac dispersion of the surface states existing 
at the interfaces. An axial vortex is included, piercing the slab along the normal to the interfaces 
with the superconductors. It is known that, when proximity is s — wave, quasiparticles in the gap 
are Majorana Bound States, localized at opposite interfaces. We report the full expression for the 
quantum field associated to the midgap neutral fermions, as derived in the two-orbital band model 
for the TI. When proximity involves odd — parity pairing, midgap modes are charged Surface An- 
dreev Bound States, and they originate from interfacial circular states of definite chirality, centered 
at the vortex singularity and decaying in the TI film with oscillations. When the chemical potential 
is moved away from midgap, extended states along the vortex axis are also allowed. Their orbital 
structure depends on the symmetry of the bulk band from where the quasiparticle level splits off. 
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I. INTRODUCTION 

Topological insulators (TIs), bulk insulators, with 
metallic surface states protected by time reversal invari- 
ance, are attracting widespread interest for potential ap- 
plications in nanoelectronics and spintronics^i^. Three- 
dimensional (3D) TIs, such as Bi2Se^ or Bi2Te^^, 
have two dimensional (2D) metallic surface states, whose 
band structure consists of an odd number of Dirac 
cones, centered at Time Reversal (TR) invariant mo- 
menta in the surface Brillouin Zone (BZ). Surface sen- 
sitive experiments such as angle-resolved photoemis- 
sion spectroscopj^ and scanning tunneling microscopy 
(STM}^ have confirmed the existence of these exotic sur- 
face metal states with a Dirac-like energy dispersion. 

Boundary states can be characterized with the help 
of minimal models for the various geometries, whose 
corresponding Hamiltonians are classified generalizing 
Altland-Zirnbauer symmetry classes^ and the topologi- 
cal invariants in mapping from the configurational space, 
parametrized by the k— vectors, to the real space. Sur- 
face Dirac fermions can become chiral when time rever- 
sal is broken by (e.g.) ferromagnetic covering and they 
can be bound to defects such as magnetic domain walls. 
The zoology of all the possible boundary states available 
has been presented in full in Ref.s lOillOl. A model which 
successfully describes e.g. the (1,1,1) surface of the (oth- 
erwise centrosymmetric) Bi2Se^ close to the T point of 
the surface BZ has been introduced in Ref. 4 and can be 
taken as the prototype model for TR invariant TIs. It is 
based on the continuum limit for long wavelengths and 
includes appropriately the spin-orbit interaction which is 
the crucial feature for the topological protection. It is a 
two band model involving two species of orbitals, of even 



and odd parity, with a bulk gap 2M. 

We consider a TI film terminated dX a,n x — y plane, 
sandwiched between two superconductors. Here we 
choose z perpendicular to the quintuple layers^i. The 
metallic states at the interfaces become superconducting 
by proximity effect. The nature of the induced super- 
conductivity when the chemical potential is within the 
bulk gap is unknown. Solving this puzzle is extremely 
important in view of the possible existence of Majorana 
Bound States (MBSs) localized at interfaces between con- 
ventional superconductors and TIs^. It has been argued 
that induced singlet, s — wave superconducting correla- 
tions, with order parameter Ag, could turn the pairing of 
the boundary states (with M > >> A^) into an effec- 
tive spinless odd — parity pairing of the p — wave type: 
Px — ipy^- In this case a vortex piercing the structure 
can bind a Majorana fermion at each interfacei^"— . It 
has been shown that, for the case of an s — wave pairing, 
in presence of an axial vortex, the zero energy Majo- 
rana Fermion can survive even when fj, is moved within 
one of the bulk bands by p- or n-doping, till the spin- 
orbit coupling is able to supply a Berry phase of tt to its 
wavefunction^^. Roughly speaking, the minus sign that 
the Majorana fermion acquires in turning once around a 
loop, which is intrinsic to its fermionic character, should 
be compensated by the spin-orbit texture, in order to 
provide a single valued solution of the eigenvalue prob- 
lem. 

Still, there is no evidence that superconductivity in- 
duced in the boundary states would be of the singlet, 
s — wave type. 

It has been found that few per cent Cu doping of 
Bi2Se3 bulk crystals moves the chemical potential /i 
within the bulk conduction band and the superconduct- 
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ing phase transition takes place at Tc ~ 3.8-K'— li^. Odd- 
parity pairing seems to be favored by strong spin-orbit 
couphng^S. Topological odd-parity superconductors are 
fully gapped in the bulk but have gapless Surface An- 
dreev Bound States (SABS). It has been noticed recently 
that a mirror symmetry (composed of the odd mirror 
reflection in the yz plane x ^ —x and of a Z2 gauge 
transformation A — — A) characterizes the helicity of the 
SABS^''^'^. Superconductivity induced in the boundary 
states of a TI slab by a singlet superconductor, odd with 
respect to the mirror plane z = Q could be of the " polar" 
type, as explained below, with an order parameter pro- 
portional to Pz and a vanishing spin projection onto the z 
spin quantization axis, which is pinned to the normal to 
the slab by the spin-orbit interaction. Strictly speaking, 
such a superconducting order would have a nodal line of 
excitations exactly at the x — y plane, which allows for 
the sign inversion at the boundary. However state of the 
art TI flakes are known to be plagued by impurities and 
vacancies in the bulk and it is likely that bulk impurity 
bands arise with little dispersion which provide a non 
vanishing Fermi momentum pp also in the z— direction. 
Therefore, we argue that the proximity gap never closes, 
even when /i is located within the bulk gap. Otherwise, 
the presence of the gap nodal line at the surface would 
affect the lifetime of the states bound at a vortex core, 
but not their topological origin. 

In this paper we explore the nature of the zero en- 
ergy excitations corresponding to quasiparticles bound 
to an axial vortex piercing a sandwich S/TI/S geometry 
and carrying a magnetic flux. While topological super- 
conductors can be well classified when TR invariance is 
preserve d^^i^^ , by the breaking of TR symmetry the topo- 
logical protection can be washed out. We compare the 
two kind of superconducting induced orders in the limit 
Af > /i >> A: a) an s-wave induced order parameter, b) 
a pz-wave (polar), odd parity order parameter. 

The two-orbital model appears to be adequate for both 
the cases. Using this model we find that a proximity su- 
perconducting gap opens at fi located in the Dirac cone 
dispersion of the boundary states. A vortex piercing the 
slab binds zero energy excitations inside the supercon- 
ducting gap depending on the value of /u. However, their 
nature is strongly dependent on the type of supercon- 
ducting ordering, a) or 6), induced by proximity: 

a) Majorana zero energy excitations are bound to the 
vortex as expectedi^ (called MBS henceforth). A 
neutral quasiparticle is localized at each of the in- 
terfaces with the superconductors and we present 
an analytical explicit expression for the quantum 
field that describes the excitation in the two-orbital 
model in the simple case of /i = (midgap state). 
Orbital parities are mixed and spin is not con- 
served. The two spacially separated modes are TR 
mates, no matter that the vortex breaks TR sym- 
metry. It has been argued that even moving the 
chemical potential inside the bulk band the vortex 
could host a MBSii. 



b) In the case of odd parity topological superconductiv- 
ity, with fi inside the bulk bands, a linear dispersion 
of SABS arises in the superconducting gap, close 
to fc = 0^°. These states are MBS. In this work we 
show that when fi is in the bulk gap, a vortex hosts 
SABS, as well. However they are no longer neutral 
fermion excitations. By breaking the TR symme- 
try of the odd-parity topological superconductor, 
the vortex makes the zero energy SABS excitations 
turn into a pair of Dirac states, loosely bound to 
the vortex core. To prove our statement, we solve 
explicitely the /i = O(midgap) case, in full analogy 
with the s — wave pairing case described in a) . The 
modes, having odd and even orbital symmetry, re- 
spectively, propagate along the vortex axis with op- 
posite chirality. At /i 7^ 0, if the chemical potential 
matches an energy split off from the bulk conduc- 
tion or valence band, zero energy excitations exist, 
with mixed type of orbitals. Expectation value of 
the spin of these excitations is, in any case, zero. 

The important novelty in our calculation is the assump- 
tion of polar order parameter for fi within the bulk gap 
in the two-orbital model. Polar order, i.e. triplet pair- 
ing with zero projection of the Cooper pair spin along 
the normal to the boundary plane, is expected to be pre- 
ferred in view of strong tendency to in-plane helical trans- 
port induced by spin-orbit coupling . We fully account 
for the two bulk band structure in the TI and for the 
third dimension z orthogonal to the plane of the slab. 
The orthogonal direction is crucial for the description of 
these bound states. The celebrated argument by Fu and 
Kanei^ is based on an effective two-dimensional system 
describing the flat boundary. It is easy to see that its 
generalization is unable to account for the proximity in 
the polar state. This is what we shortly report on, in 
closing this Introduction. 

The starting point of Fu and Kane^ is the Dirac dis- 
persion of the surface electronic states, proximized by an 
s — wave superconductor. The model can be effectively 
reduced to spinless electrons with linear energy disper- 
sion in a p — wave pairing field. Here we want to extend 
their argument to induced odd parity proximization. 

The original argument can be rephrased as follows. We 
consider just an effective Dirac-like Hamiltonian for the 
surface states at the flat boundary at z = O.The velocity 
v characterizes the linear dispersion. The proximized TI 



Hamiltonian, in the basis [ipk'f, ipki, V'lfcii ^V"- 
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with k± ^ k e^'^fc. In this basis, s ~ wave pairing takes 
A, ■ 



the form Ag = ( ^ j ^^"^ ^^^^ energy eigen- 



values are ±\/ (±vk — ^) + A^. If we disregard the two 
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bands which are far off the chemical potential, it is easy to in terms of spinless operators corresponding to spinorial 
show that an effective BdG Hamiltonian can be written wavefunctions: 
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The new effective Nambu Hamiltonian reads: 
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[vk - fi] A.-e-'"^ 
A^e*^*" - [vk - fi] 
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Indeed, e.g., the matrix element , in the basis given 

above is: 
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Eq. (121) provides an effective Hamiltonian for spinless par- 



kx 



which implies an effective 
15l. a vortex can 



tides with 

p — wave pairing. According to refll 
sustain a MBS in such a system. 

On the other hand, the odd-parity mean field cou- 
pling X^ctct' AffV tpl V'^r / + ^-^^ is described by a vec- 

tor d{k) which is odd for k —k. The generic, uni- 
form order parameter is a matrix in the spin space 
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Ap = 



. In our basis, it reads: 
If d{k) cx k itself, the ef- 



(fc) = [d(k^ ■ (s ^S2) 

dz dx idy 
X ~\~ idy dz 
fective Hamiltonian is the same as in eq.®, but 

with As replaced by \d\\ \ with d\\ = (dx, dy). In the model 
dz disappears altogether. This result doesn't look con- 
vincing, because it seems to be an artifact due to the 
restriction of the model to the two-dimensional x — y 
boundary plane. By contrast, in view of the strong spin- 
orbit locking between spin and momentum in this plane, 
it is likely that a polar "p-wave" phase is stabilized, with 
vanishing z— projection of the Cooper pair spin. This 
would imply an order parameter d = (0,0, dz) oc kzZ. 
In the bulk, such an order parameter would have a gap 
with a nodal line at the Fermi surface, exactly at kz = 0, 
because Det{Ap) (x fc^. However, this may not be the 
case close to the surface where kz is undefined. Besides, 
state of art flakes of topological insulators (prepared, for 
instance, by Chemical Vapor Deposition) are not impu- 
rity free in the bulk, and we argue that a non vanishing 
Fermi wavevector in the z~ direction is anyhow present 
in the bulk, thus providing a non vanishing superconduct- 
ing gap also at the boundary surface. Given this choice 
for the order parameter, our heuristic approach leading 
to Eq.([3]) would make, in any case, the superconducting 
order parameter disappear in H^^^ ^ what is suspicious. 

This argument shows that, when odd parity symmetry 
for the superconducting order parameter is chosen, with 



vanishing spin component of the pair in the direction 
orthogonal to the surface, a full approach which includes 
the third dimension is mandatory. 

The structure of the paper is the following. In Sec. II 
we first introduce the two-orbital model in the long wave- 
length continuum limit, close to the F point, which is of 
the Bogolubov-De Gennes type-". The model applies to 
a bulk crystal with superconducting pairing correlations. 
It will be convenient to deal with the odd — parity or- 
der parameter case first, and rephrase the model for the 
s — wave pairing. This is done in Sec. HI, where we 
derive the opening of the superconducting gap at the 
chemical potential, which can be located within one of 
the two bulk bands (thus giving rise to what is called a 
topological superconductor, Ea.([TT|)). or within the bulk 
gap. In the latter case the superconducting gap opens 
in the energy dispersion of the boundary states which is 
Dirac cone like, Eq. p^ . The neutral excitation bound 
to the vortex at the slab boundary in the case of s — wave 
proximity is derived analytically for in Appendix A and 
discussed in Sec. IV. A. In Sec. IV. B we report on the 
midgap (^ = 0) charged excitations bound at the vor- 
tex core in the odd — parity pairing. Their field operator 
and its spacial dependance is derived in Appendix B. Ex- 
tended states along the vortex core are also derived, for 
the case of /i 7^ 0. Results are summarized in Sec. V. 
From our analysis, the main conclusion that we draw is 
the observation that there are no MBS with odd — parity 
superconductive proximity and zero projection of the pair 
spin, in the two-orbital model. 



II. MODEL HAMILTONIANS 

We start from the Hamiltonian introduced by H. Zhang 
et alA to describe three dimensional layered systems like 
Bi2Se3 Sb2Te^, and Bi2Te3. They use a 4-dimensional 
basis-space which is the direct product of the spin space 
and of the orbital parity space originating from pz or- 
bitals, denoted as: {bf,t), b^t), bf,i), bM)}- Here 
g{u) denotes even (odd) parity for k ^ —k. On such a 
basis, the model Hamiltonian is (we assume the Fermi 
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velocity v to be isotropic and we put it equal to unity 
here and in the following): 

ho = -McTa + i(Ji [sidx - S2dy) + i(J-iSj,dz ■ (4) 

Here Sa and Oa {a, a S {1,2,3}) are Pauli matrices in 
the spin and orbital space, respectively. The hat, over ho 
reminds of its 4x4 matrix structure. M includes second 
order derivatives M = M{z) + Cd'^^ +cd^, parallel to the 
flat x—y boundary plane terminating the quintuple layer. 
M is half the bulk gap. Non trivial topology is guaran- 
teed by the condition M, C,c > 0, which corresponds to 
the inversion of the bulk bands. 

The Bogolubov-De Gennes (BdG) mean field Hamil- 
tonian, in the presence of a superconducting pairing re- 
quires an 8 X 8 matrix structure whose compact form 
is: 

HB,G{k) = ( ) , (5) 

where ho and A are 4X4 matrices. In the BdG rep- 
resentation, even parity singlet pairing requires that 
^s'(^) = ^s{—k), while odd-parity triplet pairing sat- 
isfies Aj(fc) = -Ap{~kp.. 

In our case, given the Hamiltonian 

H - ^iN = [ho- Pl)t^ + Hpair (6) 

where the Pauli matrices act in the Nambu space, 
the matrix structure of the off-diagonal pairing Hamilto- 
nian Hpair can acquire different forms, depending on the 
order parameter symmetry. In the even parity, s-wave, 
singlet case, we define = (V'ut^uj,) = {'^g^i'gi) = 
— (V'ttiV'Mt) ~ "(V'siV^st)' assumed to be independent of 
the orbital. It gives rise to a pairing Hamiltonian 

Hpair = ~i{AsSyT+ + h.c). 

On the other hand, it has been proposed that, when 
doped with few percent Cu, the Bi2Se^ undergoes the 
superconducting phase transition with an odd-parity or- 
der parameter—. The pairing Hamiltonian for the polar 
ordering presented in the Introduction is 

H^air = {\cfySyT++h.C.), 

where Ap = (V'«tV'si) = -(V'st^ui) = -(V'siV'nt) = 
(V'uj.V'gt) tti^ odd-parity orbital order parameter. The 
change in signs in the expectation values for Ap arises 
from the fact that the operators act on a triplet pair 
with zero spin projection along z. 

III. TI BOUNDARY STATES 

In this Section we discuss the opening of the super- 
conducting gap induced by proximity in the TI, when 
the bulk gap M is much larger than the superconduct- 
ing gap A. Depending on the doping, the chemical 



potential ^ can be within a bulk band with dispersion 
rjj: K, ±vM^+l?, or in the gap of the insulator, where 
the electronic states localized at the boundaries have a 
Dirac like dispersion. Here we choose z perpendicular to 
the quintuple layers^i, which would give helical bound- 
ary states in the surface plane. Spin orientation is ex- 
pected to be in the plane, as well. The Hamiltonian for 
the two different pairing symmetries, which we have in- 
troduced in the previous Section and the corresponding 
energy spectrum is given here below. 

Odd-parity pairing 

Let us first write down the Hamiltonian for the odd- 
parity pairing. It is convenient for our purposes to choose 
the basis 

B= [■(/'gt'V'MiiV'ui ^-^^t Vt^t' V's4,,'0s4, ^-■^/'«t '^V ■ 

(7) 

A matrix acting on this basis has 4X4 blocks cor- 
responding to components with exchanged parities 
{g,u). We address these blocks with Pauli matri- 
ces denoted by T°. Within each of these blocks, 
a Nambu particle-hole component structure arises, 

of the kind: [(V't, V'4.); (V'ji — V"!)]"^: addressed by 
Pauli matrices C". The spin space is addressed by 
Pauli matrices S"' {a, a, a = 1,2,3)). In this rep- 
resentation, the particle-hole conjugation is defined 
by S = (l X X S^) IC, where K. is the complex 
conjugation of the wavefunctions. Explicitly, the ap- 
plication of S onto the basis B gives, apart for /C: 

EB = ^^«i,-'0^t,'0«t ^ -031 ^ "09-1,- V'ut] • 

Time Reversal operator is — [1 x 1 x iS'^] JC 
(G^ = —1). This transforms the basis into 

[i^ul,-1pgt,-1pg1- \ Klpgl,-1puf,-'4'ut ^ "V'gi , 

apart for IC. 

The full 8X8 Hamiltonian in the basis B looks like: 
where 

/ tM - fi i{dx~ idy) ±Ap \ 

p _ i{dx+ idy) ±M - fi ±Ap 

±A; ±M + fi i{dx-idy) 

\ ±a; i{dx + idy) tX + m / 

(9) 

(the derivatives act on the ip fields). An unitary 
transformation, by changing the basis, B^ — > 

rewrites this Hamiltonian into the BdG form of Eq.(IS]). 
It can also be shown that another unitary transformation 
maps the Hamiltonian of Eq. ([8]) into the one of Refl20!. 

A plane wave representation can be used for a trans- 
lationally invariant material in k space. It is easy to 
check that S H{k) 2-^ = —H{—k). Time reversal trans- 
formation provides Q H{k) = H{—k) only for A* = 
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—Ap. The chirality operator, F = i s 8 is such that 
TH{k)T^^ = H{k), provided A; = -Ap. 

Bulk excitation energies, A, of the Hamiltonian of 
Eq.® are given by (/c^ ^ kf. + kl): 



(10) 



± 2^(fc2+M2)(^2 + |Ap|2) + fc2^2, 

Here we have kept only linear terms in the Hamilto- 
nian by dropping the laplacian appearing inside (i.e. 
Ai — )■ M). When the chemical potential is in the con- 
duction band {fj, > M >> |Ap|), we can neglect |Ap|2 
in the square root. Defining single particle energies 



V Af 2 + fc2 the eigenenergies take the form: 



|Ap|2. 



(11) 



While two of the resulting bands are very far from ^, the 
other two describe the opening of the superconducting 
gap at Vk^ 

It is known that at an interface between a TI and a triv- 
ial insulator, there are delocalized boundary states. Their 
energy dispersion is described by a Dirac cone e = ±fc|| , 
where fcy lies on the boundary. An interface between a 
TI and a trivial insulator with the same bulk gap, but 
with no band inversion, can be easily mimicked just by 
keeping the linear terms in the derivatives of the model 
Hamiltonian of Eq. Q and by changing the sign of M be- 
tween the two half spaces2^. The resulting wavefunctions 
of the boundary states are concentrated on the plane at 
z = and decay exponentially in the z direction normal 
to the surface, with the decay length \M\/hv. In view 
of the fact that the linearized model Hamiltonian just 
reproduces the long wavelength behavior, this boundary 
condition is good enough^^ . Analytically continuing the 



bulk eigenenergies of Ea. lfTO)) with fc2 
we get: 



— ^ 



.f2 



A = ±V(±fc||-Ai) 



|Ap|2, 



(12) 



which would describe surface superconductivity with a 
gap opening at < M. Again, two of the bands of the 
excitations are very far off the Fermi energy, while the 
other two feature the opening of the superconducting gap 
in the Dirac dispersion of the boundary states. 



as explained in Section II. A, appears as an off diagonal 
contribution in an Hamiltonian of the following structure: 



where 



HI 

id^ + {C+As + h.c.) 



id, + (C+A, + h.c.) 



(13) 



H: 



I TM - i{dx- idy) 

i {dx + idy) ±M — i-i 










±M + ^ i{dx — idy) 
i{dx + idy) + 

(14) 

The Hamiltonian in the absence of the vortex has the 
properties QH{k)e~^ = H{-k) and F7J(fc)F-i = H{k), 
only if As is real. Again, if the wavefunction decays 
as exp(— |M||z|) at the boundary, the spectrum is given 
by Ea. (jl2p . with Ap — ^ As. An unitary transformation 
which changes the basis into 



UB 



-•Pgt^ i'ul, V'gJ., -V'lt I^H'V'si,V'l^, -V'gt 



(15) 

transforms the Hamiltonian in the convenient form of 
Eq.s © and (0), with Ap ^ Ag. 



IV. QUASIPARTICLE STATES BOUND AT A 
VORTEX LINE 



As shown in the previous Section, the even-parity and 
odd-parity induced proximity give rise to the same ma- 
trix form of the model Hamiltonian, but in different 
bases. While the basis B of Eq.© is appropriate for 
the odd-parity superconducting correlations, giving rise 
to Eq.s®,®, even-parity superconducting order is de- 
scribed by the same Hamiltonian matrix when the basis 
is UB given by Eq. p3|) . We now search for zero energy 
excitations corresponding to quasiparticles bound to a 
vortex piercing an heterostructure S/TI/S in the form 
of a slab laying in the x — y plane, with boundary flat 
planes at z = 0,L. Let the axis of the vortex be along 
the z axis. It is appropriate to move to cylindrical co- 
ordinates with radial coordinate, r, measured from the 
vortex singularity and azimuthal angle around the vor- 
tex axis, 9. Outside the vortex core, the Hamiltonian 
now reads: 



Even-parity pairing 

In the chosen basis of Eq.©, the uniform order pa- 
rameter As, assumed to be independent of the orbital. 



H= 



where 



id. 



C3 id. 



(16) 
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±A* e 




±M - n 
±M + n 





-i e-^' {dr - Ide 
TM + /i 



2r; 



Here g = ±1 is the charge of the vortex and ~ Tiv/A 
is the radius of the vortex core. A stands for AgOr Ap, 
depending on the actual superconducting order. We have 
implemented sign changes with respect to Eq.sdS]),®, to 
take care of the fact that all of the derivatives should act 
to the right hand side. We have also added the vector 
potential associated to the vortex, which, far away from 
the vortex core, takes the form of a pure singular gauge: 



A. s — wave, singlet proximity 



= 0, Ae{r)^--dex 
r 



X = , (18) 



{(f> = hc/2e is the flux unit). The phase factor e'^'^^ breaks 
the TR invariance, which holds when A^ is real ( Ap is 
purely imaginary). The procedure to search for zero en- 
ergy eigenstates of the Hamiltonian of Eq. (IT6l) is sketched 
in Appendix A for s—wave pairing and in Appendix B for 
odd — parity pairing. In the next Subsections we report 
the results. 



When proximity induces s — wave, singlet supercon- 
ducting correlations, an axial vortex of charge q = ±1, 
binds Majorana quasiparticles at the interface with the 
topologically trivial superconductors. The zero energy 
eigenstates are found along the lines sketched in Ap- 
pendix A, by matching solutions inside and outside the 
vortex core. Its boundary is defined as a circle of radius 
hv/ A. An analytic derivation of the quantum field 
in the two-orbital model can be given far away from the 
vortex core in the limit oi fi — (midgap MBS). Two 
zero energy real fermion fields, localized far apart at the 
two boundary surfaces of the slab, z ^ 0+, , take the 
form, outside the vortex core: 



7(z^0+) (X e-^' K,/2{A,r) { 



-ie/2 ic 



e e"'" -(Agt 



-i0/2 iqt 



j{z ^ L ) (X e 



-A (L-z) 



K 



1/2 



(A,r) [-Z 



' - e-^o/' e-' t] + e/2 ^^ ,e ^^^^ ./2 . « t] | , 



with A |M| and z > in the TI. Here K± are the 
modified Bessel functions: i^±i/2 (wr) = e~^'^ j ^fwr , so 
that the excitations are localized also in the surface plane 
and the wavefunction is, of course, normalizable. The 
decay length scale is ^ £_o- 

The two Majorana excitations mix both u and g or- 
bitals and are not eigenstates of the spin. It is shown in 
Appendix A, Ea.s (|A6|) . that ^{z ~ 0) and 7(2 ^ L) form 
a time reversed pair, notwithstanding the fact that the 
vortex breaks TR. This is a confirmation of the fact that 
they are neutral excitations, not influenced by the pres- 
ence of the magnetic field. Having a e**^/^ factor, they 
change sign when moved along a loop about the vortex 
singularity, as expected. 

Inside the vortex core, the solution requires an r- 
dependent order parameter As[r) together with the cor- 
responding vector potential A{r) and should be matched 
with the one given previously at r = ^q. Solutions for 
any ^ within the bulk gap, require complex decay lengths 
A = Ai -|- 1X2 in the z— direction. One can impose zero 



wavefunction at z = and use the fact that, being the 
wavefunction associated with ijjua odd in z, a prefactor 
like the e"^^^ C0SA2Z can appear, while, in the case of 
ipgi;, which corresponds to an even z— function, the pref- 
actor should be of the type e~^^^ sinA2Z. 

To derive MBSs at finite < Al , all the vector com- 
ponents should be involved and the solution can be found 
only numerically. At finite /i, states extended along the 
vortex core can also be conceived. They are also ex- 
tended on the boundary surface, in the form of circular 
waves centered at the vortex singularity. In fact, in case 
the solution is of the form e^'^ ^ , it is easy to see that 

the Bessel functions solving these equations have to be 

fi 2^ 

of the first kind, i.e. J±i/2 (Agr), or H^^^'^ i'^sr), which 
describe normalizable circular states, delocalized in r. 
This is briefly shown at the end of Appendix A. These 
states, which are not neutral fermionic states, could be 
scattering states of the kind derived in Reffl^, to study 
coherent transport. Analogous states occur when prox- 
imity is odd — parity and they will be presented at length 
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in Appendix B for that case. The appearance of these 
states may be responsible for the "vortex phase transi- 
tion" , which is expected to destroy the MBS locahzed at 
each of the boundaries in z = 0, L^. 

B. Proximity induced odd-parity pairing 

An approach similar to the one of Subsection A can 
be used in the case of the Hamiltonian of Eq. p^ . with 



odd — parity pairing and zero spin projection along the 
spin quantization axis which is pinned at the normal to 
the surface of the slab. In this case, the basis is given by 
Eq. (I?]) . In a film geometry a midgap solution of the kind 
of the one given in Eq. ([TS]) , localized at the interface and 
bound to the vortex, can be found. However this is not 
a MBS. In Appendix B we report the derivation of this 
zero energy eigenstate in full analogy with the s — wave 
proximity case. One of these SABS at = is: 



4-; 



(r > lo, 9, z) cx e-"'' {i wr) ■ { 



,i(l+<2)e/2 



-j(l+g)e/2,,t 



r 



(20) 



Here z > in the TI (M, C, c > 0) and we have defined 
Ap = iA' with A' real. We get: 

— in — — fli + i 02, iw = ai — i c a2/C 

ai = ^J{C~ c){M + CA'2) - C2A'V(C - c). 



02 = 



CA' 

C-c 



M 



[c-cY 



M,C,c>0, C-oO. 



This excitation only involves fields referring to u orbitals. 
By choosing the complementary set of non vanishing 
components for the vector solution, a SABS arises, which 
involves the ipga fields only. For r > ^o, the function 
decaying in z has also an oscillatory component. The 

(2) 

function of r is a Hankel function H^^!^{i wr) of complex 
argument and has an oscillator factor exp— i oi, as well 
as a decaying exponential factor exp —caijC. For r < ^o, 



^(kV) ~ il')^'l^[K! r)-^H}^i^{K! r) is the combination of Han- 
kel functions that converges to zero at the origin (i.e. the 
point where the order parameter vanishes). In our "hard 
core" approximation, the value of is fixed by matching 
the two solutions of Eq. (pn|) at the core boundary. 

These behaviors qualify the result as a SABS, which, by 
inspection, is not a MBS, but a Dirac fermion. It is not 
an eigenstate of TR and it decays along the vortex core, 
in an oscillatory fashion. By using the projector Pi, = 
(1 - r)/2 onto the left (L) chiral state (F = lEB defined 
in Section HI), it is easy to check that the combinations 
given here in the asymptotic region out of the vortex 
core, are L— chiral at 6* = 0, i.e. of the form: [V'ui + 
i i^l^Y] and ['i/'ti-j- — i V'L ] . The partner state to the one 



(2), 



given in Ea.pop. which involves the i\jg„ field operators, 
is the _R— chiral mate combination at = 0, i.e. : [V'gf -I- 

* ^gj.]' \^gi ~ * '/'gtl-The spin expectation value for these 
quasiparticle states vanishes. However, it is interesting to 



note that just one spin component has a non vanishing 
angular momentum around the vortex line. This is, in 
Ea. (l20p . the down spin for q = 1 or the up spin for q — 
— 1, respectively. A similar feature is found in half vortex 
excitations of the ^He A— phased. 

Away from the midgap, levels can be found for Dirac 
Fermion excitations, which correspond to circular waves 
propagating at the interface inward or outward the vor- 
tex singularity and merging into the film by travelling 
across the slab, along the vortex line, with a radial local- 
ization length hv I Is! . These states involve both u and g 
orbitals. The U~^{ni — 0) ^'(r, z) field, of spacial depen- 
dence ~ e"*^ K^i{A'r), for r > £_o in the topologically 



non trivial slab, is derived in Appendix B [ Ea. (IB16p ]. 
The inverse length scale k ss \/Jl/c can be found in 
Eq.i 



(IBlSp . The location of these quasiparticle levels is 
derived by matching the solution inside the core to the 
one outside i t^'^i^^ . The matching fixes the value of fi 
at which these excitations imply no energy cost. Similar 
levels split off the conduction or the valence band and 
reside in the bulk gap. Their orbital structure depends on 
the symmetry of the bulk band from where the quasipar- 
ticle level splits off. 

We have checked other choices for the non vanishing 
vectorial components with no success and we conclude 
that there is no possibility for a MBS to exist, when 
proximity induced pairing is odd — parity. 



V. SUMMARY AND DISCUSSION 

Topological insulators hold the promise for future de- 
velopments in low power spintronics and in quantum 
computing. It has been argued that Cu doped Bi2Sez 
could become an odd-parity topological superconductor 
when /i moves within the conduction band, because of 
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dopingi^. In contrast, doping seems to be unable to in- 
duce superconductivity in Bi2Te^. Even in the case 
when the superconductors are conventional metals with 
s — wave pairing, it can be questioned whether the pair- 
ing in the TI is even — parity (singlet) or odd — parity 
(triplet) in nature. Therefore, it is interesting to char- 
acterize the nature of superconducting proximity at the 
interface S/TI—. We have considered a TI slab within 
the two-orbital model which has been successfully intro- 
duced to describe band inversion in various TI, particu- 
larly Bi2Se3 and Bi2Te^. The two bands are made of 
orbitals which are even {tpga) and odd (ipua) with respect 
to the surface plane with the normal oriented along the 
z-axis , which terminates a quintuple layer at z = 0. The 
model Hamiltonian used here is applicable both to a full 
topological superconducting state and to a TI with su- 
perconducting pairing correlations induced by proximity. 
We perform the analytical matching of the state at the 
interface between the TI and a superconductor of trivial 
topology essentially by accounting only for the the first 
order z— derivative appearing in the linearized Hamilto- 
nian. Hence the matching is performed just by impos- 
ing the continuity of the wave-function at the interface. 
In the case of an exponentially decaying wave-function 
on the topologically non trivial and topologically trivial 
sides of the interface, this is the same as changing the sign 
of the mass M at the boundary. This simplified method 
was shown to provide results^^ which are not inconsistent 
with a full treatmen1iiii2i. However, matching conditions 
at the interface z = are a minor concern in this context, 
as we have dealt with both situations, s — wave proximity 
and odd — parity superconducting pairing, on the same 
foot. In any other respect, our calculation includes sec- 
ond order derivatives, because they influence drastically 
the localization of the states in the planar dimensions 
parallel to the slab. The orbital shape of the even— parity 
order parameter is a combination of (V'st'^/'si) ^^'^ '-'^ 
{'4'u-\'>Pui) , while the odd— parity one, similar to a "polar" 
p-wave, involves the fields ^ga^'^ua' multiplied together 
and zero 2;— component of the total spin of the pair in 
the odd-parity case ("polar phase") and we expect that 
the spin-orbit coupling pins the z— spin quantization axis 
parallel to the normal to the slab surface. This feature 
could favor the proximity at the interface with a singlet 
superconductor and makes the comparison between the 
even — parity singlet and the odd — parity superconduct- 
ing pairing more intriguing, in particular in the search for 
MBSs. In both cases, superconductivity opens up a gap 
at the chemical potential, which, in our case is immersed 
in the Dirac cone dispersion of the boundary states. We 
have considered the case of an axial vortex piercing the 
slab along the normal to the surface plane and we have 
derived analytically the nature of the excitations within 
the superconducting gap, with a direct comparison of the 
two superconducting orderings. Accounting for the coor- 
dinate z in the quasiparticle wave functions, has allowed 
us to find out whether the bound state is localized at one 
of the surface boundaries, or it is an extended state along 



the vortex axis, across the slab thickness. 

Our main concern is to give an answer to the question: 
when the chemical potential is located in the bulk gap, 
can MBS appear in the two-orbital TI model, localized 
at the vortex singularity and squeezed at the interfaces 
with topologically trivial superconductors, for any type 
of proximity ordering? The answer is negative, as we 
briefly sum up, here below. 

In the case of the s — wave superconductive proximity, 
MBSs exist, localized at the vortex singularity. This fact 
is known since the work by Fu and Kane^. It is also 
known that MBS could survive if the chemical potential 
penetrates the bulk bands ( > . Here we recover 
the same result for < A/, by adopting the two-orbital 
model. There is one neutral (Majorana) fermion bound 
to the vortex, localized at each interface (top/bottom) 
with the superconductors of trivial topology. They de- 
cay exponentially away from the interface and do non 
hybridize, if the slab is thick enough. We give the full 
expression of the fields for the simplest case that can be 
handled analytically, i.e. = 0, Eg. dTSl) . They mix both 
orbitals ipga^'^ua' but they do not mix spin projections. 
It is interesting that the two partner MBS, localized at 
opposite interfaces, form a time reversed pair, notwith- 
standing the fact that TR symmetry is broken by the 
vortex flux. This fact points to the neutral nature of 
these states. 

Away from /i — 0, all the components of the vector 
spinor are involved and a simple analytical solution can- 
not be exhibited. However, by examining our analytical 
derivation it is easy to argue that, if we move the chemi- 
cal potential to non zero values, the states turn out to be 
extended Dirac fermions along the vortex axis or propa- 
gating waves ingoing or outgoing the vortex_ singularity. 
Similar solutions have been studied in Ref,16' as the lim- 
iting case of two MBS at opposite interfaces hybridizing 
signiflcantly, till they eventually delocalize. 

In the case of odd — parity superconductive proximity, 
the vortex is unable to bind neutral excitations at the 
interfaces. To prove our conclusion, we put again the 
chemical potential at the midgap, /i = 0, and we follow 
similar steps as the ones that led us to confirm the pres- 
ence of MBSs in the even — parity case. States localized 
at the two interfaces are SABS, originating from Dirac 
fermions of opposite chirality. Being SABS excitations, 
they are localized close to the boundary surfaces and de- 
cay inside the slab in an oscillatory way, in the region 
outside the vortex core (see Eg. ipD]) '). They are also lo- 
calized close to the vortex core, again in an oscillatory 
fashion. Within the vortex core, as well as in the topo- 
logically trivial superconductor, they have a simple expo- 
nential decay. The inverse localization length is oc 
while the inverse wavelength of the oscillations is ruled 
by Ap. The expectation value of their spin is again zero, 
but, in this case, vector components appear of different 
spin labels, while g and u orbitals are not mixed. This 
is the reverse of what happens in the even — parity case, 
when spin labels are separated, but different orbitals are 
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mixed. It is interesting to note that just one spin species 
in the vector has a non vanishing angular momentum 
around the vortex Une, depending on the charge of the 
vortex. Being the pair a triplet spin pair, the quasiparti- 
cle excitations can acquire features of the '^He quantum 
liquid. In Appendix B, we also derive scattering solutions 
away from midgap, in the form of circular waves localized 
at the interfaces and propagating from one interface to 
the other, by ingoing or outgoing the vortex singularity. 
Their orbital structure depends on the symmetry of the 
bulk band from where the quasiparticle level splits off. 

Having checked other possibilities, we conclude that 
there is no chance of having protected Majorana excita- 
tions with odd — parity superconductive proximity and 
zero projection of the pair spin, in the two-orbital model. 
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Appendix A: Zero energy modes for a TI slab with 
s — wave superconductivity induced by proximity 

The matrices of Ea.([T71) can be rotated by us- 
ing the projector on angular momentum eigenvectors— 



Um{q,0) 



-imO 



, to get rid of the 



explicit 9 dependence. Modes appear in pair of states of 
angular momentum ±to except for m — 0. Therefore we 
must restrict our search for Majorana modes to the sub- 
space at TO = 0. We also get U^^{q,9).j\7^jjUo{q,9)jj = 

{dr + ^) +d^, independent of j. Thus ( we put hv — 1 
for the time being). 



U^Hq,0)H±Uo{q,9)^ 

TM - /i id, 
i d, 
A' 




±M-fi 


A* 



A 



±M + 

-i dn 




with derivatives dp, dm including the vector potential, 
acting on the right hand side. Outside the core of the 
vortex, in the asymptotic region, they coincide, as they 
become: 



2r 



7r\ d„ 
2r 



dr 



(1 



^ + f •(A2) 



It appears clearly that, well away from the vortex core, 
where the choice of Ag given in Ea. (|T5)) holds, the vector 
potential drops out of the Hamiltonian. 

Let us search for a zero energy eigenvector of 
the Hamiltonian of Ea. ([T5)) of the form f{z) ■ 

Ki,?72,'?2>'7il -?7i,6,-'72>Ci] ■ where f{z) decays on 
both sides away from z — and ^, rj's are complex func- 
tions of r. If ^2, ^21 ^2, V2 ^I's taken to be zero, the eight 
equations that we get reduce to: 



-{M+ -id^rji 


= 0, 


idm^l + As T][ 


= 0, 


-idmv'i + 


= 0, 




= 0, 




= 0, 


-idm Vl + 


= 0, 


-id,n - A* f]i 


-0, 


{M id, rj[ 


0. 



(A3) 



If As is real, a solution can be obtained with f{z) 
e~^^ (2; > and A > 0). We only report the p. — 
case which is analytically straightforward: = —^1 







--^-1/2 (As r) and r\x 



Ki/2 {As r), where 



K± are the modified Bessel functions: Kj-1/2 (wr) 
e~'^'^ l^pwr. The latter functions, together with the ex- 
ponentially decaying /(z), are also eigenstates of the op- 
erator 



M 



Af + C 92 + — 9^ - — r 



2r 



4j,2 



cdl 



(A4) 



with eigenvalue A ^ M, that should be thus consistently 
determined. According to the given basis of Eg. lfTSl) . the 
linear combination of fields is a real fermion {z > 0): 



. (Al) cx e ^ ^ ii:±i/2 (A, r) • { [V'gt + V'gt ^] + i ['/'«t " V-^t 1 } 



The opposite choice: Ci = = '7i = ^1 = 0, provides 
one solution decaying at the opposite side of the slab 
z ^ L, in the form: e^^'^-'^^^K Undoing the U rotation 
and gauging away the vector potential, we obtain two 
real fermion fields. They mix both orbitals and have 
vanishing expectation value of the spin projection. The 
spinor part is of the form, outside the vortex core: 



7(z ^ 0) 



_p-i9/2 i7r/4 iq 



^gf, 0, 0, e-"/4 e-*«V„t ^ e"' e^^/^ e*«V«t- 0> 0, e'^/^ 6""/^ e-'^VsT 1 (A5) 
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We have explicitly included the field operator, from the 
basis of Ea. (fTS|) . and we have put a slash in the middle 
of the spinor to highlight the internal structure of the 
states: 7(2 ~ 0) — ^ [a|/3] and 7(2 L) — >• [a'|/3']. It can 
be seen that the a' component and the (3 one are related 
by TR, as well as the a component and the /3' one. This 
is to underline that the two Majorana states form a TR 
pair, notwithstanding the fact that the vortex breaks TR 
symmetry. By adding the field components together, we 
get the final expressions or the MBSs reported in Eq.s 
(HH). 

At finite fj, a solution of Eq.s (|A3p is also possible, prop- 
agating along the vortex axis as f{z) ~ e^*''^. Again we 
require M[S,,rj\ = to satisfy Eq.s (|X3]) 1,4,5,8. How- 
ever, at difference with the ^ = case worked out above, 
in which ^'^ = —^1 were real and 771 = —rj'^ were purely 
imaginary, all the vector components will now be real. It 
follows that the Bessel functions solving Eq.s (jA3|) 2,3,6,7 
have an argument which is purely imaginary. This shows 



that they turn into extended waves ~ J±i/2 i^sf) or 
f 1 2) 

^ ±1/2 (^s'')' ^ mentioned at the end of Section IV. A. 



Appendix B: Zero energy modes for odd — parity 
superconductive proximity 

a) SABSs at 11^0 

Here we tackle first the midgap, zero energy eigenstate 
in close correspondence with the /i = MBS of the s — 
wave pairing, presented in Appendix A. We will show 
that the solution is a Dirac fermion, that can be qualified 
as a SABS. Following the same procedure as for the s — 
wave proximity, we start from the Hamiltonian of Eq. (|16p 
in cylindrical coordinates. A generic form of the solution 
is: 



771 e 



' Jx 



(Bl) 



The equations solved by the zero energy mode are: 



-{M + fi)rii e~'^ + idp^i e'^x - e'x 



A*r]i e-'^ - {M + /i) C'l e*^' - idmV'i " idzC2p^-'^' 
A*a e-*x' -I- idp^[ e'^' - (A/ - fi) r][ e»x + id^r^'^e'^^ 

-idz£.i e-*x' „ i^„^r]2e'^ - {M + £_2e'^' - Arj'^e-'^ 
id^X'i e'""' + A*r]2e'^ - {M - /i)e^e-*x' - id^-q'^e^'^ 
idzV'i - A*i2e'^' - idpC'2e''^' + [M + fi) ry^e'^x 

I 



0, 
0, 
0, 
0, 
0, 
0, 
0, 
. 



(B2) 



We search for a real eigenvector. As Ap ~ (V^gt^fJ.) the 
choice X = x' is consistent with the odd — parity pairing, 
as it gives a purely imaginary Ap = i A' ( which defines 
A' as real). By posing f{z) = e"**^^ ( with complex 
k), the appropriate "would be MBS " requires that, e.g.. 
Ill ^ Vi — £.2 — £.2 — 0. In this case, from Ea.s (|B2|) we 
get: 



2: {M-^l)£l=0 5: -{M-^i)m=0 
3: -{M+fi)£[^0 8: (7W + /i)?/^ = 0, (B3) 



which implies that, at ;U = 0, the surviving vector com- 
ponents should be eigenstates of the operator A4 with 



zero eigenvalue. Besides: 

1 : idp^i -iA'Ci-KT]2 =0 

4 : i dp£[ ~ i A' £_i + K 7]ii =0 

6 : (9m 772 - i A' 7/2 - K ^1 =0 

7: -id,nV2~^^' V2 + kCi =0 

A possible choice is £,1 = rj2 = irj and ^ — r/, -q^ 
This implies that the equations become: 



(B4) 



dpi i (k 4- A') 7? = 0, 
ap77 - i + A') C = 0, 

a,„r/ - i (k + A') e = 

d-mi + i (k + A') ^ = 0. 



(B5) 
because. 



Rewriting d — dr + l/2r in place of dm,p, 
asymptotically out of the vortex core, the two opera- 
tors coincide, a good solution out of the vortex core. 
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satisfying convergency for large r is: 77 ~ H^^^,^{—i wr) 



and ^ 



H 



(1) 



12^— i wr) with w — [k + A'). 



/2 

The latter 



turns out to be a complex parameter. Observing that 
[d - H_ 



(1) 



_l_^^2 l"^^) — 0' constraint M. ^^±1/2 



= 



is fulfilled, provided A/ + C w 
(ai = (a + a*)/2, 02 = (a - a*)/2): 



C K 



0, what gives 



-ai + i a2, 



-fli ica-ijC 



ai = V(C-c)(M + CA'^)-C2A'7(C-c), 
CA' 



02 



-, M,C,c>0, C-oO. 



Accordingly, f{z) is a decaying function of z, but it 

also has an oscillatory component. The function of r can 

(2) 

be put in the form of the Hankel function (* ^'') ^^'^ 
has an oscillator factor exp —iai, as well as a decaying ex- 
ponential exp— ca2/C. This behavior qualifies the state 
to be a SABS, which is not a MBS, however. In fact, be- 
ing the wavefunction of the form [ 0, i f , ^, 0\i ^, 0, 0, — 
in the basis of Eq. ([T]) the field operator of the excitation 
is: 



H 



1 (« wr) ■ I 



-iTr/4 



+i 



+ e'"/^ Vut ^] } ■ (B6) 



which correspond to a mid gap SABS. 

Inside the vortex core we do not require an exponential 
decay in r, but a zero of the wavefunction at r = 0. Let us 
define f{z) = e^'^'^ for z > 0. In this case the Eq.s (|B2)) 
become: 



[tdp^i -iA'(,'i + V2] = 

[i dp^[ - z A' a - V'2] = 

-« dmri2 - i A' 772 + i k' 6] =0 

a„i77^ - i A' ?72 - i k' ^J] = 0. 



(B7) 



Here dp^m do not have the form of Ea. (jA2p . because, the 
vector potential does not take the asymptotic expression 
typical of a singular gauge within the vortex core. Also 
A'(r) is expected to have a linear r— dependence close 
to the vortex axis24. A simple way to approximate the 
eigenfunction is to consider a hard core with A' = and 
to overlook the vector potential difference. By putting 
?/2 = ^1 = ^ and ^[ ~ ri2 = rj we get the two equations: 

l(-7): d^ + K'T^ = 
4(-6): dri-K'^ = 



l2 

'K rj. 



(B8) 



Appropriate solutions of Eqs. ijBSp arc the Hankel func 

(1 2) 

tions H^^^2{K'r). Besides 



M K, r?] = -^M-{C- c) k' 



(B9) 



what defines k' as a real parameter for M, C — c > 0. 
Having two possible normalizable solutions, we can im- 
pose that their combination vanishes at r = 0. We thus 
obtain ^(r) and r](r): 



eM = <i/2K) + <i/2K), 

riir) = -i{H%2i»r)-H%2i<^'r)} 



(BIO) 



which are real for any r. They converge to zero at the 
origin and are such that rjir) — ^(r). In conclusion, inside 
the core, the vector is: 

^e-«'^ [0,e(r),-r;(r),0|-77M,0,0,^(r)]^. (Bll) 

In the basis of Eq.Q, this gives: 

- e-^'^ {{4>ui- i^ui - (V«t + V^ut ^) } -(612) 

Undoing the transformation J7(m = 0) we get Eq. ipO)) of 
the text. Our "hard core" approximation depends on the 
parameter ^o, which fixes the core boundary. The value 
of can be determined by matching the inside and the 
outside solutions at the core boundary. 

By choosing ?7i , jyi , ■^2 , £,'2 ^lo^i zsi'O ^.nd the other spino- 
rial components zero, the partner state of Ea. pop . in- 
volving the fields, can be found. 

b) extended states along the vortex line 

Here we show that a Fermi Dirac state solution is 
also possible, involving real vector components, which 
describe a wave travelling along the vortex , with 
z— dependence e"*^ and k real. This is a bound state 
within the bulk gap. All the ^'s, rj' s have to be non van- 
ishing. If we pose ri[ = -772, 773 = and = ^2, ^2 = 
^1, the Eq. (jB2p can be solved with: 



(B13) 



6 - i^i/2(A'r), 6 = A'-i/2(AV) 
771 = ifi/2(AV), 772 = if_i/2(A'r) 

by requiring the following : 

1: + ^)r7i + K7/2 = 5: - /x)?72 - k ?7i = 

provided 



1': ap6-A'^; = o 



5': 9p6-A'e^ = 



The first two are compatible, if, after substitution of the 
eigenvalue to the operator A^, the determinant : M.^ — 
-\- = Q vanishes. Similarly: 

8: -(A1 + Ai)?72 + K7?i =0 4: -(X - ^)77; - k 77^ = 
provided 



8': dp£'2-A'£2 = 

provided 

6': ~id,nm + ^' V'2^^ 
: -{M + ^l)£,[+Ki'2 = Q 
provided 
3': 9™r;; A' 771 = 



4': dp£[-A'£i=0 
2: {M- fi)£i + K 6 - 

2': a™77i + A' r;l = 
7: iM^^l)C2 + ^iCl^0 

7': dmV2+^'V2 = 
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The value of k for any /i is fixed by solving the determi- 
nantal condition: 



+ {hv nf - =Q (B14) 



( fiv has been restored). By posing m = M + CA' we 
get: 



K = 




m {hv) (hv) 



4c2 



(B15) 



which gives a real k, because^^: in > (fiv)^ /2c. The 
plus sign in front of the square root of Eg. ljBlsp has been 
chosen on physical grounds, in order to obtain that the 
wavelength of the propagation along z increases when 
/i goes deeper down in the bulk gap. Qualitatively is 
At ~ [ij c. The matching between the inside and the 
outside of the vortex core will fix the value of /i at which 
the excitation is zero energy^S. Hence the location of the 
energy level is fully determined. 



Just outside the slab, at the interface with the topo- 
logically trivial material is m < 0, so that < and 
the solution decays with z away from the surface. Being 
K purely imaginary, the vector components of the eigen- 
function acquires alternatively an i factor, as is the case 
leading to Ea. (jB6p . This changes the Bessel functions 
K±i/2 into the corresponding Hankel functions. While 
in the case of = these functions were localized in r 
anyhow, because of complex argument, in this case they 
are delocalized in r, because the argument, A'r, is real. 
This implies that, at the interface, the waves become cir- 
cularly propagating inward or outward the vortex line. 



The consequence is that this state is a traveling wave 
along the vortex axis in the non trivial topological ma- 
terial , while it decays outside at the interface with the 
trivial material and propagates outward or inward at the 
surface. Its energy is localized in the gap. This is not a 
Majorana state, however. Using Eq. ljBip for the vector, 
we obtain: 



which, in the basis of Eq.([7]), provides: 

^ U-\m = 0) *(r, z) cc e'^' if±i (A'r) • [(V^^t + ^at ^) + {^H - i^m ^) + (^«T + + {^9l + ^)] -(216) 
The final step would be undoing the transformation U{m — 0^ 
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